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Abstract 

Recent progress in the fabrication of quantum dots using silicon opens the prospect of observing 
the Kondo effect associated with the vaUey degree of freedom. We compute the dot density of 
states using an Anderson model with infinite Coulomb interaction U, whose structure mimics the 
nonlinear conductance through a dot. The density of states is obtained as a function of temperature 
and applied magnetic field in the Kondo regime using an equation-of-motion approach. We show 
that there is a very complex peak structure near the Fermi energy, with several signatures that 
distinguish this spin- valley Kondo effect from the usual spin Kondo effect seen in GaAs dots. We 
also show that the valley index is generally not conserved when electrons tunnel into a silicon 
dot, though the extent of this non-conservation is expected to be sample-dependent. We identify 
features of the conductance that should enable experimenters to understand the interplay of Zeeman 
splitting and valley splitting, as well as the dependence of tunneling on the valley degree of freedom. 
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I. INTRODUCTION 



Experimentation on gated quantum dots (QDs) has generally used GaAs as the starting 
material, due to its relative ease of fabrication. However, the spin properties of dots are 
becoming increasingly important, largely because of possible applications to quantum in- 
formation and quantum computing. Since the spin relaxation times in GaAs are relatively 
short. Si dots, with much longer relaxation times , are of great scientific and 

technological interest 5|. Recent work has shown that few-electron laterally gated dots can 
be made using Si 6| 7|], and single-electron QDs are certainly not far away. 

There is one important qualitative difference in the energy level structures of Si and 
GaAs: the conduction band minimum is two-fold degenerate in the strained Si used for 
QDs as compared to the non-degenerate minimum in GaAs. Thus each orbital level has a 
fourfold degeneracy including spin. This additional multiplicity is referred to as the valley 
degeneracy. For applications, this degeneracy poses challenges - it must be understood and 
controlled. From a pure scientific viewpoint, it provides opportunities - the breaking of the 
degeneracy is still poorly understood. 

Previous experimental work has shown that in zero applied magnetic field the splitting 
of the degeneracy is about 1.5 ± 0.6 /leV in two-dimensional electron gases (2DEGs) and 
that it increases when the electrons are further confined in the plane [8]. Surprisingly, the 
splitting increases linearly with applied field. Theoretical understanding of these results has 
historically been rather poor, with theoretical values much larger than experimental ones for 
the zero-field splittings |9|[l0| and theory also predicting a nonlinear field dependence 11|. 
^■lecent work indicates that consideration of surface roughness may resolve these discrepancies 
. In Si QDs, a related issue is also of importance: what is the effect of valley degeneracy 
on the coupling of the leads to the dots? The spin index is usually assumed to be conserved 
in tunneling. Is the same true of the valley index? 

Recent measurements of the valley splitting in a quantum point contact show a valley 
splitting much larger than in the 2DEGs, about 1 meV in quantum point contacts [s^. A 
QD in a similar potential well is expected to produce a valley splitting of the same order of 
magnitude. The overall picture of the degeneracy is as shown in Fig. [1] for a single orbital 
level of a Si QD. Of particular interest is the fact that a level crossing must occur, and the 
rough value of the applied field at this point is Bcr ~ 2.5 T, given a reasonable zero-field 
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FIG. 1: Schematic diagram of the energy levels of a single orbital in a Si QD in a magnetic field. 
The quantum numbers are even (e), odd (o) valley states and up and down spin projections for 
spin 1/2 states. A level crossing occurs at the magnetic field Bcr- 

valley splitting of 0.5 meV, and the valley-spin slope of magnetic field for the "o t " state 
to be 0.1 meV/T, as it is in Hall bars. 

The aim of this paper is to study the Kondo effect in the context of transport through a 
Si dot. The Kondo effect was originally discovered in dilute magnetic alloys. At low temper- 
atures, the electron in a single impurity forms a spin singlet with electrons in the conduction 
band, thus causing an increasing resistance as the temperature is reduced to zero. Since the 
first observations of the Kondo effect in GaAs QDs, there has been considerable experimen- 
tal and theoretical work done, mainly because QDs provide an excellent playground where 
one can tune physical parameters such as the difference between energy levels in QDs and 
the Fermi level, the coupling to the leads, and the applied voltage difference between the 
leads. The level of scientific understanding of this purely spin Kondo effect in QDs is on 
the whole quite satisfactory. The basic phenomena are as follows. At temperature T < Tk, 
where the Kondo effect appears, a zero-bias peak in the dot conductance is observed. An 
applied magnetic field splits the peak into two peaks separated by twice the Zeeman energy 
2gfiBB. These linear peak energy dependencies on the applied magnetic field have been 
observed in GaAs QDs [l^. The Kondo effect only occurs when the occupation of the dot 
is odd. Clearly, Si dots will have a much richer phenomenology. There are multiple field 
dependencies as seen in Fig. [H and the additional degeneracy can give rise to several Kondo 
peaks even in zero field. Furthermore, the Kondo temperature is enhanced when the degen- 
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eracy is increased. We build on previous work on the orbital Kondo effect. For example, an 
enhanced Kondo effect has been observed due to extra orbital degree of freedom in carbon 
nanotubes (CNTs) 1^ 16] and in a vertical QD, with magnetic field-induced orbital degen- 



eracy for an odd number of electrons and spin 1/2 [17|. More complex peak structure in 
their differential conductance suggests entangled interplay between spin and orbital degrees 
of freedom. We also note that a zero-bias peak has been observed in Si MOSFET structures 

Our goal will be to elucidate the characteristic structures in the conductance, particularly 
their physical origin, and their temperature and field dependencies. Once this is done, one 
can hope to use the Kondo effect to understand some of the interesting physics of Si QDs, 
particularly the dependence of tunneling matrix elements on the valley degree of freedom. 

This paper is organized as follows. In Sec. II, we determine whether valley index conser- 
vation is to expected; in Sec. Ill, we introduce our formalism; in Sec. IV, we discuss the 
effect of valley index non- conservation; in Sec. V, we present our main results; in Sec. VI 
we summarize the implications for experiments. 



II. TUNNELING WITH VALLEY DEGENERACY 



The basic issue at hand is the conductance of a dot with valley degeneracy, and the 
influence of the Kondo effect on this process. The dots we have in mind are separated 
electrostatically from 2DEG leads on either side. These 2DEGs have the same degeneracy 
structure. Thus before we tackle the issue of the conductance we need the answer to 
a preliminary question: is the valley index conserved during the tunneling process from 
leads to dots? Non-conservation of the valley index will introduce additional terms in the 
many-body Hamiltonian that describes the dots, and, as we shall see below, it changes the 
results for the conductance. In this section, we investigate the question of valley index 
conservation in a microscopic model of the leads and dots. We do this in the context 
of a single-particle problem that should be sufficient to understand the parameters of the 
many-body Hamiltonian that will be introduced in Sec. III. 

If we have a system consisting of leads and a QD in a 2DEG of strained Si, then both the 
leads and the dot will have even and odd valley states if the interfaces defining the 2DEG 
are completely smooth. The valley state degeneracy is split by a small energy. If the valley 
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index is conserved during tunneling, then the even valley state in the lead will only tunnel 
into the even valley state in the dot. On the other hand, if the valley index is not conserved, 
then the even valley state in the lead will tunnel into both the even and odd valley states 
in the dot; similarly for the odd valley state in the lead. We can estimate the strength of 
the couplings between the various states in the leads and the dot by calculating the hopping 
matrix elements between them. 

Consider a system of a lead and a dot in a SiGe/Si/SiGe quantum well, separated by a 
barrier of height T4. If H = 00, there will be no tunneling between the lead and the dot 
and the eigenstates of the whole system can be divided into lead states |^^(e,o)) dot 
states l^^(eo)) where e, o arc the even and the odd valley states. Call the Hamiltonian for 
this case Hq. When Vb is lowered to a finite value, there will be some amount of tunneling 
of the lead wavefunction into the dot. Call the Hamiltonian for this case H. This can be 
thought of as a perturbation problem where the n*'* perturbed eigenstate \n) is no longer 
the n*'* unperturbed eigenstate \-nP) but acquires components along the other unperturbed 
eigenstates The perturbation Hamiltonian is H' — H — Hq. H' considered as a function 
of position is not small, but the matrix elements of H' are small. This means that in the 

order perturbation theory we can write the perturbed n^^ state as 

N = K) + ElO^^ (1) 

where T4„ = (A;°|if'|n°) is the hopping or tunneling matrix term. 

If we further assume that the tunneling Hamiltonain conserves spin and that only one 
orbital state need be considered on the QD, we can expand the perturbed even valley lead 
wavefunction (with tunneling) as 

Usually, in a perturbation problem, the unperturbed wavefunctions and the hopping matrix 
elements are known and used to find the perturbed wavefunction. However in the present 
case we use a model to compute the perturbed and unpertubed wavefunctions and use these 
to determine the unknown hopping matrix elements Ve,e Ve,o- These will be needed in later 

sections. 

From the above equation, 

Ve,e = (|i^L-i^Del)(*^el^Le) (3) 
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Ve,o = (|i?L-^?)ol)«l*Le) (4) 

The perturbed and the unperturbed wavefunctions and the corresponding energies are cal- 
culated by using an empirical 2D tight binding model with nearest (vz) and next nearest 
neighbor interactions (uz) along the z (growth) direction and the nearest neighbor interac- 
tions (vx) along the x direction (in the plane of the 2DEG). Potential barriers and edges of 
the quantum well are modeled by adjusting the onsite parameter (e) on the atoms. This 
is an extension of the 2-band ID tight binding model outlined by Boykin et al. jo], lo| . 
considering only the lowest conduction band of Si. This is a single particle calculation with 
the magnetic field 5 = 0. Note that the parameters from the Boykin et al. model are chosen 
precisely so as to reproduce the two-valley structure of strained Si. 
The single-particle Hamiltonian for the system is 

H = Y1 [e\(i){x,z)){(P{x,z)\+v^\(P{x,z)){(P{x + l,z)\+Vz\(i){x,z)){<i){x,z + l)\ 

X 

+ Uz\(l){x,z)){(j){x,z + 2)\ +fx.|0(x- l,z)){(j){x,z)\ +Vz\(p{x,z- l)){(f){x,z)\ 

+ Uz\(P{x,z-2)){<j){x,z)\ ] (5) 

The parameters defining the system are as follows: the lead is 80 atoms long, the barrier 
between the lead and the dot is 16 atoms wide along x and 0.3 eV high, the QD is 38 atoms 
(^ 10 nm) wide along x and 38 atoms wide along z, the barrier along z defining the quantum 
well is 20 atoms wide on each side and its height is 0.3 eV. The nearest and next nearest 
neighbor interaction terms along z are Vz = 0.68264 eV, and Uz = 0.611705 eV, and that 
along X is Vx = —10.91 eV. The schematic of the system is shown in Fig. [2] 

We consider two kinds of 2DEGs, one with no interface roughness at the interface and 
the other with a miscut of about 2°. We find that the valley index is conserved in the case 
of smooth interfaces. We model the miscut interface as a series of regular steps of one atom 
thickness and width defined by the tilt, e.g., a 2° tilt corresponds to a step length of about 
30 Si atoms. In this case, the valley index is no longer conserved. Strengths of the coupling 
between the e — e and e — o valley states is found to depend strongly on the relative distance 
between the edge of the closest step and the edge of the quantum dot. In Fig. [31 we plot the 
hopping matrix elements as a function of this relative distance. The terms show oscillations 
that have a period equal to the step size in this case. 

The coupling strengths are also found to depend on the tilt of the substrate on which 
the 2DEG is grown. We show in Fig. HJ the hopping matrix elements as a function of the 
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FIG. 2: The system of the lead and the dot in a miscut quantum well with the potential profile 
along the x-direction. The open circles represent Si atoms sandwiched between the SiGe atoms in 
the barrier region represented by the filled circles. 

stepsize (oc 1/ tilt) and find that the hopping terms vary rapidly at small step sizes {i.e large 
tilts) but slowly at larger stepsizes [i.e., smaller tilts). 

In a quantum well, the phases of the rapid oscillations of the electron wavefunctions (along 
the z-direction) are locked to the interface and in a well with smooth interfaces, the phase 
remains the same for a particular valley throughout [12]. But in a miscut well, the phases 
are different for electrons localized at different steps even for the same valley eigenstate. 
The resultant valley splitting is due to interference of phase contributions from all the steps. 
Hence the valley splitting can be very different in the lead and the dot. The same argument 
explains the conservation of the valley index across the barrier in the lead-dot system. In the 
case of a smooth interface, the e eigenstate in the lead has same phase of the z-oscillations as 
the e eigenstate in the dot, and is 90° out of phase with the o eigenstate in the dot. Hence, 
there is a finite overlap between the e — e wavefunctions, while the e — o overlap cancels out 
due to interference. This preserves the valley index during the tunneling. 

In the case of a miscut well, we can no longer strictly label the eigenstates e or o, but let 
us continue to do so just for convenience. Here, because e in the lead does not really have 
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distance between step edge and Qdot edge(atoms) 



FIG. 3: (a)Tunneling matrix elements (K,e,K,o) in units of eV and (b) Ratio of the tunneling 
matrix elements (Fe,o/K,e) as a function of the relative distance between the step edge and the 
quantum dot edge for constant barrier height and width (for a miscut QW). The solid line with 
triangular markers indicate the Vg^e term and the dashed line with open circles indicate the T4,o 
term. Here the stepsize is constant at 30 atoms 2° tilt) 

a single phase along z but is composed of different phases at each step and similarly for o 
in the dot, there is a fair chance that there will be a non-cancellation of the phases between 
the e, o wavefunctions depending up the the extent of tunneling. This will give rise to a 
finite hopping matrix term between all the dot and lead states and the valley index o, e can 
longer be said to be conserved. 

The results show that the conservation or non-conservation of the valley index depends 
sensitively on fine details of the dot morphology such the proximity of a step edge to the edge 
of the dot on atomic scales. It is not likely that good control at this level can be achieved, 
and thus some sample dependence must be expected. This means that it will normally be 
necessary to include hopping terms that do not conserve valley index into the Hamiltonian 
in order to understand the Kondo effect in realistic Si QDs, where the interfaces are rarely 
smooth. 



8 




FIG. 4: (a)Tunneling matrix elements (V"e,e,K,o) in units of eV and (b) Ratio of the tunneling 
matrix elements {Ve,o/Ve,e) as a function of the stepsize for constant barrier height and width (for 
a miscut QW). The solid line with triangular markers indicate the Ve,e term and the dashed line 
with open circles indicate the T4,o term. Here the step edge coincides with the quantum dot edge. 

III. EQUATION-OF-MOTION APPROACH IN THE U ^ oo LIMIT 

Sec. II has given us some insight into the nature of the tunneling between the leads and 
the dot. We now introduce a Hamiltonian to describe the full many-body problem, and 
discuss the computational method we will use to solve it. This Hamiltonian must include 
the single-particle energy levels of the leads and the dot, the tunneling matrix elements that 
connect these levels, and the Coulomb interaction between electrons on the dot. We shall 
use the Anderson impurity model: 

ikma "^"^ ikma 

ikma m'u'^ma 

Here a indicate two-fold spin 1/2 and m two-fold valley indices, m means the opposite 
valley state of m. The operator cli,^^{cikma) creates (annihilates) an electron with an energy 
£fe in the i lead, i e L,R, while the operator f^f^ifma) creates (annihilates) an electron 
with an energy Ema in the QD, connected to the leads by Hamiltonian couplings Vo,ife and 
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Vx,ik which correspond to an electron tunnehng between the same valley states and different 
valley states, respectively. We assume that the Vo(x),ifc do not depend on spin index a. U 
is the Coulomb interaction on the dot. It is assumed to be independent of m. Previous 



studies on the orbital Kondo effect have mostly assumed conservation of orbital indices 



20| 



4) to SU{2) Kondo 
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19| or in double 



Recent theoretical work has found that a transition occurs from SU 
effect by allowing violation of conservation of orbital index in CNTs 
QDs 2l|. For Si QDs, we have shown from the previous section that tunnelings between 
different valley states have to be taken into account. 

Experiments measure the current / as a function of source-drain voltage Vsd- The differ- 
ential conductance G = dl/dVsd is given by the generalized Landauer formula 3| 

I = r / dw[fL{w) - fR[w))- —lm\grna{w)\ 

where D is the bandwidth and Ti^rna{u)) = 2n J2k(yoik + ik)^i'^ ~ ^k)- The f's are Fermi 
functions calculated with chemical potentials /i^ and ^r with fiL = fJ'R + ^^sd- In the fol- 
lowing argument, we assume a flat unperturbed density of states and Vo,ifc = Vq, Vx,ik = Vx- 
We then define F = TL,maiw) + Tr^^„{w) = nV^/D with = + With these ap- 

proximations we may differentiate this equation and we find dl/dVsd ~ — Ini[^m(T(eVsrf)]/7r = 
density of states(DOS). Here Q^^(w) is the retarded Green's function: 



poo 

« fma, fL »= / e^"* < {fmAt), fLiO)} > dt. (6) 







where a = w + i6. Here {, } and <, > denote the anti-commutator and statistical average of 
operators, respectively. Hence our task is to compute Qmaiw). Although the approximation 
of a fiat density of states is not likely to be completely valid, we still expect that sharp 
structures in Q^„(w) will be reflected in the voltage dependence of dl/dYgd- 

Several approximate solution methods for this type of Hamiltonian have been used suc- 
cessfully: numerical renormalization group, non-crossing approximation(NCA), scaling the- 
ory, and equation-of-motion (EOM) approach. In this paper, we use the EOM approach to 
investigate the transport through a Si QD at low temperature, in the presence of magnetic 
field. The EOM approach has several merits. The most important for our purposes is that 
it can produce the Green's function at finite temperature, and works both for infinite-L'^ 
and finite-[/. We basically follow the spirit of the paper by Czycholl [231] which gives a 
thorough analysis of EOM in the large- iV limit (A^ is the number of energy level degeneracy) 
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and provides a suitable approximation for finite-t/ systems. Though in the present case is 
4, the Green's function obtained from EOM should give a reasonable picture in the case of 
a weak hybridization and relatively large U. The approximation starts from the equations 
of motion for Q.^^[w) in the frequency domain: 

w<^A,B >=< {A, 5} > + < [A, ff], 5 >=< {A, 5} > + < A, [H, B] > . 

This can be obtained by integrating Eq. [6] by parts. The detailed derivation for the Green's 
function Qma is given in the Appendix. We truncate higher-order Green's functions by using 



a scheme adopted by |24j . l25l . |26| . The higher-order Green's functions are of the general 
forms <^ f~^cc, ^ and -C /c"'"c, ^ (we drop indices for the sake of generalization) 
which we decouple by replacing them with 

< /+CC, /+ > = < f+c X c, /+ > 
« /c+c, /+ » = < c+c >« /, /+ » . 

We combine Eq. lA.ll IA.2t IA.171 and IA.18] and take the limit U — > oo. We obtain coupled 
Green's functions expressed by 

R{z)gm. = I- <^l> - H Ai{z) + P{z)Mfna 

R{z)Mrna = < fLUa > +F ma{w) + P{z)g^., (7) 

where M.fha =^ /mo-, fma ^ the correlation function between different valley states. We 
denote z = w — Ema + ^ij^ma and z = w — Ema + £p^fha- The propagators R{z), R{z) and 
P(z), P{z) are defined as follows 



R{z) 
















R{z) 






- E '?'(^) 








p^ma 


P{z) 


= SfeH(i 


- E A(2)) - Srf(w) - ima{w) + 'La{w)Fjna{w) 












P{-z) 












p-^ifia 





Note that / 7^ mcx, p 7^ ma. The various functions in R{z), R{z), P{z), and P{z) are defined 
in Eq. IA.19tlAl28l in the Appendix. Moreover, they can be self-consistently computed. We 
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proceed to solve for Qma and M.ma- 



Qr, 



R{z) - P{z)P{z)/R{z) 



^ ^ < fmafrna > +F^^{w) + (/ - El^ma < rij > - J2l^^a Ai{z))P{z) / R{z) 

R{z)-P{z)P{z)/R{z) ■ ^' 

P{z), P{z) express the correlations between "ma" and "ma" states. If Vx = 0, then 
P{z) = 0, Coupled Eqs. Hare decoupled and reduced to the Green's function with conserved 
valley index. Furthermore, when Vx 0, peak structure will alter since the function T,d{w) 
in P{z), P{z) has a second-order logarithmic divergence at the Fermi level. Now we have a 
set of self-consistent equations for Q^u from Eq. [HI since Ama{z), < rima >, < fmafma > and 
other terms such as Ima{w), Fma{w), Jma{.z) are integrals over Q^a and A^mo- themselves. 
We solve by taking an initial guess for and M.mcr of the forms 

1 - T.l^ma <ni> 



MiJw) 



W - Ema + iNV-K 
^ fmafma > 



These have the expected structure of a broad peak of width A^F (A^ = 4) centered around 
the energy Ema for Qrna and Emu for Aima, and the narrow Kondo peak(s) around the Fermi 
level. We then iterate to self-consistency, at each stage determining the occupation numbers 



< > and expectation values < f^afrna > from [271] 

< rima > = [ dwfFD{w)lmgrnc 

TT J 

< fmafrha > = / c/w/fD (w)Im7W, 

TT J 



TT 

Here fpn is the Fermi function. 



IV. EFFECTS OF NON-CONSERVATION OF VALLEY INDEX 

Consider all four-fold valley and spin levels to be degenerate, meaning no magnetic field 
and no zero-field valley splitting. We have R{z) = R{z), P{z) = P{z) and z = z = w. We 
define the Kondo temperature as the temperature at which the real part of the denominator 
of the Green's function vanishes. As discussed by Hewson [28], this definition gives the 
temperature that controls the thermodynamics. We rewrite Eq. [8] in order to explicitly 
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express two contributions from Qma, 

^ ^ 1 - El^ma <ni> - E;^m^ < f^^Ua > +Fma{w) 

^""^ 2{R{w) - P{w)) 

1 - E;^ma <ni> - El^ma M^) " (< f^gfrna > +^ma(w)) 

2{R{w) + P{w)) 

We hence find two Kondo temperatures corresponding to the solutions of 

ReR{w) + ReP{w) = 0, and 
ReR{w) - ReP{w) = 0. 

To second order in the V^s, these equations are 



(10) 



w 



£m.a -2^fFD{£k)— ^ = and 



ik 



W -£k 
2 



W 



0. 



ik ^-^k 

In general each equation has three solutions, but only one of them lies above the cut —D< 

w < on the real axis. To find T^, we set w = ksTx, Ep = and assume ksTx <^ D. As 
the temperature tends to zero, T^'s satisfy 



knT, 



B-I-Kl 



TT 



kBTK2 = Sma In 

TT 



D 

kBTK2 



D 



where 



TT 

'd 

D 



2V^^ + 2V^ + {Vo + Vxf 
2V^ + 2V^ + (K) - Vxf 



Define Vq — V cos 0, Vx — V sin (p, p — sin 20, with 0</3<l. P — implies valley index 
conservation, while P — 1 means that tunnelings between same valley states and different 
valley states both take place. 
We express Fi, r2 as 



r2 = 



D 
ttV^ 
D 



(3 + /?) = (3 + /3)r 

(3-/3) = (3-/3)r 
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(a) ^ = (b) /? = 0.5 (c) /3 = 1 

FIG. 5: We plot the DOS with various /3 values. Other parameters are A = 0, T = 0.2 meV, 
D = W meV, and Ed = —2 meV. In each graph, the red dot-dashed curve is plotted from the 
first term of Eq. [TOl corresponding to Tk2, and the blue dashed curve from the second term, 
corresponding to Tki- The purple solid curve is a combination of the two. (a) (3 = 0, T/Tki{P = 
0) = 0.6, and TxiifS = 0) = 2.83 x 10"^ meV. The red dot-dashed curve overlaps with the 
blue dashed curve as Tki = Tk2 with the peaks located at Tki{P = 0). (b)[(c)] (3 = 0.5 [1], 
T/Tki{(3 = 0.5[1]) = 0.6. The blue dashed curve shifts right away from the Fermi level and its 
peak sits at TxiiP = 0.5[1]) = 1.26 x 10"3[3.88 x 10"^] meV, while the red dot-dashed curve shifts 
left towards the Fermi level with its peak located at Tx2(/3 = 0.5[1]) = 3.49 x 10"^ [1.51 x 10^*^] meV. 
The inset in (c) shows the peak corresponding to Tk2{P = 1)- Note that the blue dashed curve is 
more pronounced than the red dot-dashed curve because Tki > Tk2 as /3 > 0. 

with r = nV"^ / D. Assuming fc^^i^i, kBTK2 -C Emu, we obtain 

ksTKi = D exp (7r£^./(3 + P)T) (11) 
kBTK2 = D exp (7r£^./(3 - /3)r) (12) 

There are two important observations from the above result: first, the Kondo tempera- 
ture depends exponentially on /3; secondly, there are two Kondo temperatures which split 
when /3 > 0. For the case /? = (conservation of valley index) we find Tki = Tk2 = 
Dexp (Tre^o-ZSr), as expected from the EOM approach. When < (3 < 1, meaning either 
Vo > Vx OT Vo < Vx, electrons can now tunnel through same valley states or different valley 
states and the two Kondo temperatures Txi(/3) and Tk2{P) split. When /3 = 1, we have the 
maximum Tki{(3 = 1) = (DT^2(/3 = 1))'^' > Tk{(3 = 0) > Tk2{P = 1). 

This enhancement of the Kondo peak and hence the Kondo temperature is caused by 
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interference between two tunnelings, one with valley index conserved, and the other without. 
This is evident from the fact that the T^-'s depend on the relative sign of Vq and Vx- 
Constructive interference increases the Kondo temperature whereas destructive interference 
decreases it. This /5-dependent Kondo temperature is somewhat different from that in fiol . 
where the authors use slave-boson mean field theory. It would be interesting to investigate 
the difference between these two approaches. It should be mentioned that Tk{P = 0) 
acquired by the EOM approach underestimates the true Kondo temperature which is 

kBT^ = Dexp (7r£:„^/4r) 

This is certainly a deficiency of the EOM approach, and it has been discussed in the literature 



2J]. The underestimation of the Kondo temperature from the EOM approach leads to its 
failure in determining, for example, the local spin susceptibility in the Kondo regime j26| . 
However, as we see in the following, the disappearance of Kondo peak structure suppressed 
by the temperature and the magnetic field is only scaled by the Kondo temperature Tki(2){P) 
and hence the EOM approach faithfully captures the Kondo effect. 



V. DENSITY OF STATES 



In Fig. 1 we showed the level structure of a single orbital on the QD. The four energy 
states suffer both Zeeman splitting and valley splitting. We now formally define these ener- 
gies as Ema = £d + {'^/'^ + fJ'vB){6m,e - Sm,o) + gfJ'BB{6„^^ - 6„^l), whcrc Ed is the bare energy of 
the dot and B is the applied magnetic field. A is the zero-field valley splitting. Experiments 
in 2-DEGs also show that the valley splitting increases linearly with the magnetic field Q]. 
For a QD the valley splitting slope depends on the size of the dot. A small Si QD can have 
fly relatively comparable to gfiB ~ 0.1 meV/T, with g = 2. 

In the following, we consider four different situations with different values of the three 
parameters B, /3, A and in each discuss the peak structure in the DOS in full detail. We 
have chosen the parameters F = 0.2 meV, D = 10 meV and Ed = —2 meV. This leads to 
rather low Kondo temperatures, but is favorable for illustration purposes. 
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A. B = and A = 0. 



This is the case considered in the previous section. A = is probably not achievable 
in experiments, but this case is still important to understand, since it allows us to isolate 
the effect of valley index non-conservation. We showed in the previous section that there 
are two /3-dependent Kondo peaks coming from the two terms of Eq. [TOl the first term 
corresponding to Tk2, the second term corresponding to T^i- The DOS can be similarly 
partitioned, and we plot each piece separately together with their sum. Fig. [5] shows the 
Kondo peak structure at /5 = 0, 0.5, 1. We plot each graph at temperature T = 0.6 Tki{13). 
When /? = (valley index conserved), T^i = Tk2, the two Kondo peaks overlap with the 
peak position at T^i, each contributing to half of the total weight. As (3 increases and 
tunnelings between different valley states become allowed, Tki{I3) increases with /5, while 
Tk2{.P) decreases, according to Eq. [11] and [121 Since we set the temperature T slightly smaller 
than Tki{i3), but much larger than Tk2{P = 0.5, 1), we expect to see one pronounced peak 
shift left and sit at energy Tki{(3), and the other much suppressed peak shifts right towards 
Tk2{i3), as shown in Fig. [3 Therefore, when tunnelings between different valley states are 
allowed, it will enhance exponentially one Kondo peak with the Kondo temperature T^i, 
and suppress exponentially the other Kondo peak with the Kondo temperature Tk2- Tki 
overshadows Tx2 except in the neighborhood oi (3 = 0.5, where two peaks are relatively 
well separated and both sizable. As stated above, it is probably difficult to achieve A = 
experimentally. However, even if A is small, we expect an enhancement of the upper Kondo 
temperature, making the features in the DOS associated with this easier to observe. In 
this sense. Si is more favorable than GaAs for observation of the Kondo effect. 

B. B = 0, /? = 0. 

We now consider what happens when there is a finite valley splitting at zero field. We 
find that there are three peaks in the DOS (and hence in the nonlinear conductance) instead 
of the usual two. There is a central peak at zero energy and side peaks at ±A. This is 
shown in Fig. [6] The schematic diagram on the top left shows the conduction processes 
corresponding to the central peak and and the left-side peak. As clearly demonstrated 
in the schematic, there are basically two energetically distinguishable types of transitions, 
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FIG. 6: We plot the DOS with /3 = 0, 5 = and A = 0(red sohd curve), A =6 TxiiUue 
dot-dashed curve), and A =12 Tft:i(brown dotted curve) at a temperature T = 0.3 Tki where 
Tki = 2.83 X 10~^ meV. The other parameters are the same as in Fig. [5j Note that at A = 0, the 
peak shifts shghtly away from the Fermi level and locates itself at Txi- As A increases, the peak 
splits into three peaks at the positions —A, 0, A. The central peak comes from spin exchange at 
the same valley state (no energy penalty) and shifts back to the Fermi level, as expected from the 
conventional spin-1/2 Kondo effect. The other two side peaks result from inter- valley processes 
(energy cost A). A schematic on the top left depicts two transitions that contribute to the central 
peak and left-side peak. 

the inter-valley transition and the spin-flip transition at the same valley. Because valley 
degeneracy is already broken with a splitting A, the inter-valley transitions cost an energy 
penalty A, therefore shifting the peaks away from the Fermi level to ±A. They represent the 
transitions where an electron in the conduction bands enters into one valley state while an 
electron hops out of the QD from the other valley state (red dotted curve in the schematic). 
Furthermore, there is no energy cost in the spin- flip transition at the same valley, since 
spin degeneracy is still preserved, so the corresponding peak locates itself at the Fermi level 
(blue solid curve in the schematic). This is basically the same scenario as traditional spin- 
1/2 GaAs QDs. Fig. [7] shows that increase of temperature suppresses these three Kondo 
peaks equally, demonstrating that they are scaled by the same Kondo temperature. 

We further calculate the Kondo temperature with A 7^ 0. Since /3 = 0, we reduce the 
real part of the denominator of Qma to 

W-Ema- J2 ^c,/(^) = (13) 
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FIG. 7: We plot the DOS with (3 = 0, A = 12 Tki at various temperatures: T/Tki = 0.3 (blue 
solid curve) , T/Tki =0.6 (red dot-dashed curve), T/Tki =1.2 (green dashed curve), and T/Tki = 
1.8 (brown dotted curve). The Kondo temperature T^i = 2.84 x 10~^ meV. The other parameters 
are the same as in Fig. [5l Here the three Kondo peaks are equally suppressed by increase of the 
temperature. 



Consider "ma" = "o j " energy state with the energy Eq,^ =6^ — A/2. AsT^O, Eq. [13] 
can be rewritten as 



w - Ed + A/2 + rin 



+ 2rin 



w + A 



D 







Again we set w = ksTxiA). In the region ksTxiA = 0) ^ | A| ^ D, because Tk{A 7^ 0) < 
Tk{A = 0), kBTxiA 7^ 0) ^ A. Therefore, we neglect linear w and the last logarithmic 



term 2rin 



D 



2rin 



. We obtain 



ksTKiA) 



Dexp{{Ed- A/2)/T -2ln\A/D\) 
{kBTK{A = 0)f 
A2 



(14) 



Thus using EOM, we find Tk{A) ~ l/A^. Eto in his paper 29|] evaluates theoretically 
the Kondo temperature in the QDs with two orbitals and spin-1/2 as a function of A, the 
energy difference between the two orbitals. He considers the case (3 = so the two Kondo 
temperatures coincide. Using "poor man's" scaling method, he finds that within the region 
where Tk{A = 0) ^ |A| ^ D, Tk{A) decreases as A increases, according to a power law 



kBTK{A) = kBTK{A = 0) ■ {kBTK{A = 0)/|A|)' 
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FIG. 8: The DOS is plotted at /3 = with the magnetic field B = (red curve), 0.00284 T 
(green curve) and 0.00568 T (blue curve). The Kondo temperature Tki = 2.84 x 10~^ meV. We 
set T/Tki = 0.2. The other parameters are the same as in Fig. [5j We set fi^ = 0.1 meV/T, 
slightly smaller than gfis = 0.114 meV/T. The central peak splits into two peaks because the spin 
degeneracy is broken, while three energetically distinguishable inter-valley transitions due to the 
presence of magnetic field split the side peak into three peaks. The transitions that result in peaks 
a, b, c, d(green curve) are shown in the schematic diagrams a, b, c, d, respectively. Moreover, each 
peak evolves with a different slope with respect to the magnetic field due to cancellation or addition 
of spin and valley slopes. The slopes of peaks a, b, c, d are 0.228, — 0.028, 0.2, 0.428 meV/T, 
respectively. 



where 7 = 1. The discrepancy of Eq. [14] from Eto's power law is again due to the deficiency 
of EOM approach to find the true Kondo temperature. 

Interestingly, as A becomes larger, the side peaks shift further away from the Fermi level 
and we recover the conventional spin-1/2 Kondo effect, with one Kondo peak at the Fermi 
level. This is in agreement with Eto's analysis 29 1. 



C. p = o. 

This is now the more general case of finite field and valley splitting, but still assuming that 
valley index is conserved. Now the peak structure in the DOS becomes more complicated. 
As already seen in Fig. [6l at i? = and A 7^ 0, there are three peaks at w = 0, ± A. By 
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applying a magnetic field, hence breaking spin degeneracy, as shown in Fig. [HI the central 
peak splits into two peaks with peak positions at ±2gfiBB by Zeeman effect. In Fig. |8l 
the schematic diagram a demonstrates a spin-flip transition that contributes to the peak 
at position 2gfiBB. The side peak above the Fermi level splits into three peaks with peak 
positions at A + 2{^y—g^B)B, A + 2fiyB, A + 2{^y+g^B)B which correspond to inter-valley 
transitions in the schematic diagrams b, c, d, respectively, whereas the side peak below the 
Fermi level also splits into three peaks with peak positions at —A — 2(/it, -|- gfiB)B, — A — 
2fi^B, — A — 2(yU^ —gfiB)B by the same token. It is easy to see that twice the Zeeman energy 
2g^BB corresponds to the energy cost for spin-flip transitions, while valley splitting energy 
A + 2fi^B corresponds to that for inter-valley transitions. The positions of these eight Kondo 
peaks have different dependencies on the magnetic field, with their slopes being expressed 
as various linear combinations of A, /i^,, and gfJ^B- It is worth noting that equally weighted 
peaks at positions ±A ±2{^^ + g^B)B, ± A ± 2(yU^ — gfiB)B are twice as small as the other 
peaks because only half the number of conduction processes contribute to the peaks. 

D. /3 = 1. 

This is the case where tunneling between same valley states and different valley states are 
equally strong. When the magnetic field is zero, there are three peaks near the Fermi level 
which we have explained in subsection B. However, as can be seen in Fig. [9], the central peak 
is more pronounced than the other two side peaks. This is due to extra contribution from 
tunnelings between different valley states. The effect of non-conservation of valley index 
can be seen more clearly when we apply a magnetic field. When B ^ and assuming the 
zero-field valley splitting A 7^ 0, we saw eight peaks since spin and valley states are no longer 
degenerate. Surprisingly, as shown in Fig. [9l when the valley index is no longer conserved, 
we not only see the previous eight peaks as in Fig. [8] when the valley index is conserved, 
but also see an extra peak at the Fermi level that is independent of the magnetic field, 
meaning it can only be suppressed by increasing the temperature. This peak comes from a 
logarithmic divergence in the function ^^(w) (cf. Eq. IA.22I) and is thus more pronounced as 
(3 becomes larger. It also indicates that through cross- valley couplings Vx, "ma" and "ma" 
form a bound state through exchanging electrons with the leads. The schematic in Fig. [9] 
demonstrates that this kind of transition through cross-valley couplings can occur without 
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FIG. 9: The Kondo peaks are shown with B = T (red curve) and B = 0.0388 T (blue curve) 
when the vahey index is not conserved. We consider /? = 1, A = 8 Tki meV and the Kondo 
temperature T^i = 3.88 x 10^'^ meV. We set T/Tki = 0.2. The other parameters are the same 
as in Fig. O Note that there remains a peak at the Fermi level even when the magnetic field is 
not zero, in contrast to the peak structure when valley index is conserved (cf. Fig. [8|). This peak 
comes from a transition where a conduction electron from an even valley state hops onto the dot 
at the odd valley state while the electron at the odd valley state hops out of the dot to a lead at 
an odd valley state, as shown in the top right schematic. This kind of transition costs no energy 
penalty. 

any energy cost. Consequently its peak sits at the Fermi level and has no dependence on the 
magnetic field. This is in some sense the most characteristic signature that could be sought 
of the valley Kondo effect, since it results from a bound state of different valley states rather 
than different spin states. 



VI. CONCLUSIONS AND IMPLICATIONS FOR EXPERIMENTS 

Silicon quantum dots show a much richer phenomenology than GaAs dots due to the 
additional valley degree of freedom. The initially fourfold degenerate level is split at zero 
applied magnetic field by interface effects and there is additional splitting at finite fields. 
We have computed the density of states on the dot in the large-f/ limit using an equation of 
motion method that has been shown to give reasonable results in related problems. This 
function should closely mimic the non-linear conductance through the dot as a function of 
voltage. 
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The nonlinear conductance near zero bias of a Si QD with a single electron is sharply 
different from that of a GaAs QD in the Kondo regime: 

1. At S = we expect three peaks near zero bias at eV — 0, ±A, where A is the valley 
splitting. The central peak is the spin Kondo peak while the other peaks correspond to 
inter-valley conduction processes. 

2. For B > 0, the side peaks further split into three peaks, the two split-off peaks corre- 
sponding to spin-flip processes and the middle corresponding to spin conserving processes. 

3. For B > 0, the central peak splits into two peaks if the valley index is conserved and 
three if it is not conserved, the latter corresponding to a valley-Kondo ground state. 

4. The temperature dependence can be complex, since there are two Kondo temperatures 
in the problem. 

5. The conservation or non-conservation of valley index depends quite sensitively on the 
precise details of the junctions and is therefore likely to be sample-dependent and difficult 
to control. 

We hope that some or all of these phenomena will be observed in future experiments. 
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APPENDIX: EOM FOR INFINITE U 

The basic equations of motion for the Green's functions are 

{W - Ema) < fma, fma'> = ^ + ^0,ife < Qfemcr, /ma > + Vx,ik < Qfema, fla > 

ik ik 

+ U J2 ^ fl^flfma, fL > ■ 

(W - £k) < Cikma, f^aa > = Vo,ik < fma, f^a > +^X,ifc < fma, f^a > • 

(W - Srna) < frna, fma^ = Yl ^0,ik < Qfcma, /ma > + Z ^X,ik < Qfema, fma > 

ik ik 

+ U E ^ fpfpfrna,fL^ ■ 
p^ma 
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{W - Ek) < Cikfha, fma > = Vo,ik < ffha, fma > +Vx,tk < fma, fma > • 

Note that l,p are shorthand for valley plus spin indices, unlike m that means only the valley 
index. Furthermore if / = ma, I = rha. The bar means exchanging only the valley index. 
From now on we simplify the calculation by changing Vo(x),ifc — * Vo(x)- From the above four 
coupled equations of motion we can obtain 

(t^ - - E ^^^) « /^<. » = i + E^^«/^-/^.» 

+ U E ^f^flfmajL^ (A.l) 
l^ma 

{W -Sma-Y. — -) < /m<x > = E < /™-' fL > 

+ E (A.2) 

The higher-order terms give rise to 

{W - Sma -U) <^ nifma, fma^ = < > + E ^o{< Uldkma, fma^ + '^ fl^Ciklfma, fma > 

^iC-ikma ) fma'> 

ik 

+ ^ ft'^iklfma, fma ^ " ^ C%,iflfma, fma ^} 

+ t/ E ^nmjfmaJma^ (A.3) 

(W - - f/) < np/^^, /+^ > = -Sp^ma < fmaffha > + E ^o{< "-pQfc^a, > 

~^ fp Cikpffhaj fma ^ ~ ^ '^ikpfpffria, fma ^} 
f^pCikmaj fma ~^ ^ /p ^ikpfmai fma 

j^p^ma 

Next, we proceed to obtain equations of motion for the higher-order functions appearing 
on the right-hand side of Eq. IA.3I and Eq. IA.4I 

{W - Ek) < riiCikma, fma > = -^l,m.a < fma^tkma > +K){< ^fma, fma > 

- E ^ '4k'lfl'^ikma, fma ^ + E ^ f l' Ci' k' iCikma ■, fma ^} 
i'k' i'k' 

+ Vx{<. riifma, /ma > - E ^ 4k'Tfl^ikma, fma > 



i'k' 

+ E f^^^i'k 'iCikmay fL »} (A.5) 
i'k' 
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(w 



^iklflfmcFi fmcr ^ 



— Sk) ^ niCikrha, fmcr ^ " ~^l,mcr < fmcr'^ikma > +Vo{<^ Ulfmcr, fmcr ^ 

- XI ^ (^k'Jl^^ikfha, fmcr > + 51 ^ fl^ Ci> k'lCikma , fmcr >} 
i'k> i'k' 

+ Vx{< riifma, /ma > - H < 4k'lfl^ikmcr, fma > 

i'k' 

+ X/ ^ fl^'^i'k'l'^ikma, fL »} (A.6) 

^fe) < f^Ciklfma, fma^ = < f^Cikl > f^ flfmcr, /ma > " ^ 4 k' iCikl f mcr , fmcr > 

i'fc' 

- E < CiH//+Q.fe,ma, /+a >} + /i+A/ma, /ma > 
i'fc' 

- E 4k'fiklfma, /ma > " E ^ CiM f,^ Ci> k'rha , fmcr >} (^-7) 

— ^/c) ^ fl^Ciklfma, fmcr ^ = < fl^Cikl > +Vo{^ ft flfmcr, fmcr ^ E ^ 4 k'l^iklfmcr , fma ^ 

- E < Wi'^Ci'fc'ma, /+a >} + Vx{^ f^flfrncr, fL > 
i'W 

~ E 4k'f^iklfma, fma ^ ^ E ^ Ciklfl^'^i'k'mcr, fmcr ^} (^-8) 
i'fc' i'k' 

l/iw + Ek-ei-ema-U) (A.9) 

l/(w; + £fc-£,-ema-f/) (A.IO) 

U < n,nj/^,, /+^ > = t/(< riiTij >+■■ ■)/{w - e^a - 2U) (A.ll) 

— < ftciklfma, fma ^ = H{< fl^ flfrha, fma ^ ~ E ^ C^k'l^^iklfrha, fma ^ 

i'fc' 

- E < CiM/;+Ci.fc,ma, fma >} + {< /i+A/ma, f'L > 
j'fc' 

- E ^ 4k'fiklfrna, fma ^ -J2 Cikl f,^ Ci' k'ma , fma >} (A- 12) 
i'k' i'k' 

— £k) ^ fl^Ciklfrha, fma ^ = Vo{^ fl^ flfrha, fma ^ ~ E ^ 4k'l^iklfrha, fma ^ 

i'k' 

- E < Cifc^/ + Q.fe,ma, /+a >} + {< /i+/Jma, /^a > 
i'k' 

- E ^ 4k'fiklf^'r, /ma > - E ^ Cikfftci'k'ma, fma >} (A.13) 
j'fc' i'k' 

l/iw + Ek-ei-ema-U) (A.14) 
l/(w + £fc -e, -ema - f/) (A.15) 

U < nzn^/ma, fL'> = U (5i,ma < n^fmafrna > + " " O/l^^ " ^ma " 2f/). (A.16) 

Note z = w — Ema + £^«7^ma and z = w — Sma + ^«7^ma- If we take f/ — > CX3, Eq. IA.9tlA.ll1 and 
Eq. IA.14tlXrT6] vanish, as do correlators And the rest involving two annihilation operators 
on the dot like <^ fi^ fifma, fma ^ also vanish as f/ — oo. After we truncate the higher- 
order Green's functions, we will encounter related integrals that by simple manipulation we 



4klflfmai fma ^ 
'^jfc;"/'/ma) fma ^ 
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transform into integrals over Green's functions <^ fma, fma ^ (Gma) and ^ fma, fma ^ 
{Airhcr)- Thus the set of EOM terminates after truncation and the Green's function Qma can 
be solved for. The equations are 



E 

ik,i'k' 



E 

ikA'k' 



Z — Sk 



Z — Ek 



" de r de' [ dw'fpoiw')^-^^''"'"''"'''''''^'' 



AnD^ J-D J-D 



z — e 



—ITT 



{z) + VoVx{B 



(z) + 



dw — + ^^(Vq B„ 



2D J-D w' -z-iS 

■Tncr^mai^Z^ -\- -Bnio-^mcr (^) ) ) • 



(z) + ViB, 



E 

ik 

E 

ik 



E 

ik 

E 

ik 



^ '^ikmafmcr > 



ma, ma 



^ ^ikmafma > 



2D 



(z)) + VxB 

ma, ma 

-2 ~ 

ma, ma 

(z)) + VxBma,ma{z)). 



Z -Sk 



with 



Ba,p{w) 



D 



z)) + VxBma,ma 



dw'fpuiw'j 

D W' — W — to 



Note that we have considered both the correlations in the conduction bands, and on the dot. 
However, conduction electrons do not mix different valley states except through cross-valley 
couplings (Vx) with the dot. Thus terms like J2i'k', 



ik 



< ^i'fc'i^''''> g^j^g order Vi3v^, while terms 



W-Ek 



'O(X)' 



such as J2ik 



W-Ek 



are of order Vq^x) ■ Eq- IA.3I and IA.4I can be thus expressed as 



{W - Ema - U - T.a{w) - T.a{z)) < Uifma, fma >= ^b{w) < Uiffna, fma ^l^ma 
+ {5l,ma{Ima{w) - T.a{w)Fma{w) + T.d{w)) + T.b{w) Ai{z)} Mfha 

+ i:a{w)Ai{z) - - Ji{z)}g^,+ < ni > +A(2;)(A.17) 

{W - Eraa - U - - < Upffna, f^a >= ^fcl^) < '^p/ma, /m<^ >p5^ma 

^p,ma ^ fma fma ^ ^p.ma^ ma 

(w) + {6p^ma{Imaiw) - T,a{w)Fma{w) + 

+Sf,(w;)y4p(z)}^m<^ + {-5p,m<^Sb(w)F„^(w;) + Sa(w;)y4p(^) - T.c,p{z) - Jp{z)}Mrna{A.18) 



where the various functions are defined by 



ik 



W-Ek 



(A.19) 
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ImAw) 



AAz) 



ImAw) 



FmAw) 



E 



2VoV: 



X 



+ Vx) < 4k'iCiki > 



(A.20) 



E 



i'k'ik 

_ 2 fD 

~2D J-D 



Z — Ek 



X 



w' — z — i6 
+VoVxiB,j{z) + BjAz))). 

ik,i'k' ^ 

2 2VoVx 
dw fFD{w ) — 

in' 



+ 



D2 



{V^BiAz) + V],Bjj{z) 



(A.21) 



+ 



l^^o B 



2D j~D w' — w — iS 

+ VoVxiBma,fhaiw) + Bfha,maiw))) . 

^0 < '^i'k'ma'^ikma > < C-i' k' ma^ikma > 



(A.22) 



E 

i'k'ik 

—in 
~D2" 



w - efc 

((Kf + V^)Bma^^Aw) + 2V^V^B^ 



[w) 



Fma{w) = 



+ iVo + V^)VoVxiB^a,mAw) + B^^^rnAw))). 
< fmaCikma > +Vx < f^aCikma > 



ik 

-2 



W - Sk 



M^) = E 



j^iiVo + V^)B„,^,^Aw) + 2VoVxB^, 

VqVx{< 4k>fiki > + < 4k'iCiki >) 



H). 



(A.23) 



(A.24) 



i'k'ik 

—iix 



Z- Sk 

r2 , T/2 



VoVxm + Vi){BiAz) + Bij{z)) + 2\/oVx(Sm W + Bii{z))). (A.25) 



E 



< f^Ciki > +Vo < f^dki > 



Z-Sk 



2D 



{{V^ + V]t)B,Az) + 2V,VxB-,Az)). 

^ 4'k'mcr'^ikmcr > < cf'k'ma^ikma > 



(A.26) 



w - 



(W) + 2Vo'V^Bma,ma 



E 

i'k'ik 

— ITT 

~D^ 

+ {Vo+V^)VoVx{B H)). 

^ ^0 < fmcr'^ikfha > +Vx < fma^ikma > 
ik 



W 



W - Sk 



2D 



{{V^ + V]i)Brn.,mAw) + 2VoVxB„ 



H). 



(A.27) 



(A.28) 



Combining Eq. IA.lt IA.2t IA.17t lA.lSt and taking the limit t/ — > oo, we obtain Eq. [7] in the 
main text. 
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